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Abstract 

Charged current events are studied at HERA with polarized beams, using a Monte 
Carlo event generator and taking into account detector smearing effects. The observed 
asymmetries are used to study the extraction of the polarized structure function g§ and 
the polarized parton distributions Am and Ad in the proton. 



■ 1 Introduction 

r- 1 . 

o ■ The flavour decomposition of the polarized quark distributions in the proton is important in- 
i 
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formation for our understanding of the spin structure of the proton. At present fixed target 
(D ". experiments the asymmetries measured in semi-inclusive data ||T| have been exploited for this 
purpose, a method which can also be used at HERA with polarized beams 0. This option 
for HERA has been proposed in ||, and leads to collisions of 27.5 GeV polarized electrons on 
820 GeV polarized protons. The high center of mass energy of HERA (y/s = 300 GeV) allows 
in addition to study charged current (CC) events which probe different combinations of spin- 
dependent quark distribution functions compared to neutral current interactions. In particular, 
CC interactions distinguish quark from anti-quark flavours. 

This paper is a continuation of the analytical calculations reported in |4[], and is based on 
a Monte Carlo study of CC events. For this purpose the Monte Carlo event generator PEPSI 
6.5 |5|, U was upgraded to contain the full electroweak structure at tree level for lepton-proton 
scattering. It was used to generate event samples which were tracked through a detector program 
to study detector smearing effects. PEPSI is a version of the program LEPTO 0, extended for 
polarized beams and targets. It contains leading order QCD matrix elements, supplemented with 
parton showers to approximate higher order effects. The hadronization in PEPSI is based on 
the Lund string fragmentation model ||. The cross section results of PEPSI have been checked 
with the analytical calculations reported in ||, and were found to be in agreement. 

The sensitivity of the CC cross section to the polarized quark distributions is shown in the 
following. The hadronic scattering tensor for lepton-nucleon interactions can be decomposed in 
three unpolarized and five polarized structure functions. The latter (gi to #5) are not indepen- 
dent. Assuming g± ~ 0, which is strictly true in the parton model, the differential cross section 
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for CC reactions on a longitudinally polarized target reads 



CPa GpMyy 



dxdQ 2 Air (Q 2 + M\ 



, ([aF^Q^-XbFs&Q^ + i^Xbg^Q^+ags&Q^S), (1) 
w) 



with x being Bjorken-x, and Q 2 the four-momentum transfer squared in the deep inelastic 
scattering process. Gf is the Fermi constant, and My/ is the H^-Boson mass; A is — 1(+1) for 
the e~ (e + ); 5 denotes the longitudinal polarization of the nucleon to be anti-parallel (+1) or 
parallel (—1) to the polarization of the lepton; a = 2(y 2 — 2y + 2) and b — y(2 — y), with y the 
usual inelasticity variable. The asymmetry is defined by 



da^+da^ aF^± hF WT 



(2) 



where in leading order 



gf = Au + Ac + Ad + As, gf = Ad + As + Au + Ac, 
gf~ = Au + Ac - Ad - As, g™ + = Ad + As - Am - Ac, 



= u + c + d + s, F? =d + s + u + c, 

F™~ =2{u + c-d-s), F™ + =2(d + s-u-c) 



(3) 



Here Ag denotes the polarized parton distribution which is the difference between the probability 
to find a parton in the longitudinally polarized proton whose spin is aligned minus the one whose 
spin is anti-aligned, and q is the unpolarized parton distribution. 

To evaluate the measurability and quality of the spin asymmetries for CC events several 
systematic studies have been performed. A total luminosity of 500 pb" 1 and a total polariza- 
tion of 0.5 (0.7 for each beam) was assumed and generated with PEPSI, compatible with the 
expectations for HERA ||. The polarized parton distributions GS-A [[| and the unpolarized 
parton distributions GRV [JT/Oj] have been used throughout this analysis. The available luminos- 
ity was distributed in equal parts for incoming electron and positron beams, and for parallel or 
antiparallel polarizations; i.e., for each case 125 pb _1 was used. However, since the proton has 
roughly two times more up than down quarks, it could be better to assign two thirds of the 
luminosity to the W + case to obtain results with equal statistical significance. Since this does 
not interfere with other polarized measurements such a scenario should be considered for future 
measurements and studies. 



2 Detector Effects and Simulation 

To investigate the measurability of the spin asymmetries a fast simulation of a HERA type of 
detector jTT], [T^] has been implemented [Tj| (for more details see for example [|14|| ). The experi- 
mental characteristic of charged current interactions at HERA is a high p t neutrino in the final 
state leading to large missing transverse momentum for the detected event. The correspond- 
ing selection of the events relies essentially on the calorimetric measurement of the hadronic 
final state. For this study we based the simulation on the calorimeter of the HI experiment. 
The simulation parameters have been tuned to real data, and some details will be given in the 
following. 

A calorimeter with a separate backward and central/forward part was considered. Backward 
here means the electron direction. The calorimeter covers the polar angle 9 (defined with respect 
to direction of the incoming proton) from 4° to 177°. The central/forward and backward parts 
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are separated at 152°. An optional acceptance in the range 10° < 9 < 170° was implemented 
to take into account the planned luminosity upgrade of HERA which would reduce the 
acceptance as a result of machine magnets which are installed close to the interaction point. 
The complete range in azimuthal angle <fi is covered. 

The influence of the inactive material in front of the calorimeter, and of losses due to the 
calorimeter edges and in the region between the central/forward and backward part was modeled 
with a set of efficiency factors which depend on the polar angle of the simulated particle. 

The smearing of the energy and reconstructed angle of the detected particles was applied as 
follows. The energy resolution of the backward part is taken to be a^ ac k/E = 0.56, and of the 
central part a cent /E = OAQ/y/E © 0.73/ E © 0.026, where the energy is given in GeV. For the 
angular resolutions <jq = 50 and cr^ = 90 (in mrad) was used. 

Further, a Gaussian smearing (with RMS 11 cm) of the interaction vertex position along the 
beam line, consistent with current HERA operation, was included in the simulation. 

The efficiency to trigger and select charged current events was set to 85±3% in accordance 
with current estimates |OJ for charged current events. The background of photoproduction 
events, cosmic rays and production of W or Z° Bosons is very small in the present analyses at 
HERA and was ignored in the following [16|, . 

We will study the asymmetries as a function of Bjorken-a;. The bin size in x as used in the 
following was chosen taking into account statistics and event migration due to detector effects. 
For each bin we have a purity]] and an efficiency^] bigger than 80%. 



3 Results 
3.1 Kinematics 

The HI collaboration selects CC events requiring a minimal missing transverse momentum, 
Prmiss-, of the hadronic final state to be larger than 25 GeV and presents results with a minimal 
Q 2 of 625 GeV 2 ||16|| . The ZEUS collaboration on the other hand has recently published an 



analysis of CC events with the selection criteria of Prmiss > 11 GeV and Q 2 > 200 GeV 2 fL7| 



The lower cut allows to extent the reach in kinematics from log 10 (x) ~ —1.5 to log 10 (x) ~ —2.0 . 

For the studies presented in this contribution Prmiss > 15 GeV and Q 2 > 225 GeV 2 have 
been used, which seem a reasonable assumption based on the present day experience at HERA. 
The results for the asymmetries, including detector effects, are shown on the left side of figure [I]. 
The error bars indicate the statistical precision of the measurement. The asymmetries are very 
large, as noticed before in H, so that the data allow for a significant measurement. The solid 
line is the result of the analytical calculation of the asymmetry using the formulae of section 1, 
and the same parton distributions as used for the PEPSI generated events. It shows that the 
detector smeared asymmetries are in good agreement with the true ones. 

Note that in our kinematic domain we have a » b, thus the asymmetry to a good approxi- 
mation directly measures the structure function ratio g^/Fi. 

w ± 

* W± « |W (4) 
1 number of events generated and simulated in a bin divided by the number of simulated events in the same 

bin 

2 number of events generated and simulated in a bin divided by the number of generated events in the same 

bin 
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The structure functions will be already precisely measured by the time the analysis of 

polarized data from HERA can be performed. This makes it possible to extract the structure 
functions from the measured asymmetries. Here we multiply the asymmetry with as 
calculated from the GRV distributions. The results are presented on the right side of Fig. [I] 
and compared with the asymmetries from the analytical calculation for g$. It shows that this 
approximation works well (to the 10-20% level) in our kinematic range. 




-2.5 -2 -1.5 -1 -0.5 -2.5 -2 -1.5 -1 -0.5 

og, (x Bi ) log, (x B; ) 



Figure 1: Spin asymmetries A w (full symbols, left side) and A w+ (open symbols, left side) 
for CC events are presented for a total luminosity of 500 pb _1 . Also shown are the structure 
functions g§ (right side) extracted from the asymmetries. The parton densities GS-A from || 
were used. The error bars represent the statistical uncertainty of the measurement. 



3.2 Energy scale of the calorimeter 

The selection and measurement of CC events depends strongly on the calorimeters of the detector. 
This implies that one of the dominant sources of systematic uncertainty of the measurement is 
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the precision on the knowledge of the energy scale of the calorimeters. The electromagnetic 
energy scale will be rather well determined for the complete detector, to better than 1%, but 
the uncertainty on the hadronic energy scale will be larger. 

To study this effect the energy scale of the calorimeters was changed in the simulation program 
by ± 10% in the backward part and by ± 4% in the central/forward region. These values are 
in agreement with those currently used by the HI experiment and thus represent a somewhat 
pessimistic scenario. We find that the effect of varying the energy scales on the asymmetries 
and the structure function g 5 is within the statistical fluctuations of the measurement, shown in 
Fig. |2|. Hence, it turns out that the effect of the energy scale uncertainty largely cancels in the 
measured ratio of cross sections (see eqn. 2). 
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Figure 2: Effect on A w± (right) and gY (left) of the variation of the absolute energy scale of 
the central/forward and backward calorimeters. 
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3.3 Acceptance of the detector 



Currently the geometrical acceptance of the central detectors at the HERA collider is about 4° 
— 177° in polar angle 9. This will be reduced by the high luminosity upgrade of HERA. For 
this upgrade it is planned to install magnets within the detectors, close to the interaction point, 
tightly surrounding the beam pipe WB. Apart from a reduced angular acceptance, the large 



amount of material added close to the interaction point will result in an area in the calorimeter 
which will be strongly affected by secondary interactions and backscattering. This area will 
have to be excluded from the analysis for CC event studies. We took a fiducial region of the 
calorimeter, which amounts to approximately 10° to 170° in polar angle, considered to be usable 
for the CC analysis. To study the effect of the reduced acceptance, the simulation program was 
changed accordingly (but the possible remaining influence of backscattering due to the magnets 
inside the detectors has not been included). The results show however that the quality of the 
measurement of A w is not significantly affected by the reduced acceptance, and is compatible 
within statistics with the result of Fig.[IJ, using the full acceptance. 

3.4 Influence of the parton density functions 

The studies presented so far have used the parton densities of GS set A. Here results for GS set C 
and the parton densities of reference JOJ for the standard version in leading order approximation 
(labelled GRSV-S-LO) are shown in Fig. |3| for A w and gY ■ ft appears that the predictions 
are very similar. This is due to the limited flavour separation possibilities in extracting these 
parton densities (from QCD fits to g\ measurements), and in the overall normalization imposed 
by the g\ measurements. However, only future measurements will show whether the present 
assumptions are in agreement with reality. 

To extract g§ (right side of the figure) the GRV parton densities were used in all three 
cases to calculate 

4 Extraction of Au and Ad 

As mentioned in section [I], we have in leading order 

gf~ = Au + Ac -Ad- As 
gf + = Ad + As - Au - Ac 



In all three parametrizations of the polarized parton densities used in this analysis Ac has 
been neglected. Also the following assumptions are made: Am = As and Ad = As. With 
these assumptions the extraction of gY + is equivalent to the extraction of Ad. Furthermore, 
gY is close to Am. This is demonstrated in Fig. f| where the ratio of Au to gY for the three 
parametrizations used in this analysis is shown. It is seen that for the kinematic range accessible 
to this measurement (log 10 (a;) > —2) all three Am distributions account for more than 85% of 
gY ■ Hence the extraction of Au can be done with a precision comparable to the gY extraction 
shown in the preceding sections. 
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Figure 3: Asymmetries A w± (right) and structure functions (left) for three different parton 
densities. In all three cases GS-A was used to calculate , needed to extract from A 1 ^. 
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Figure 4: Ratio of the Au polarized parton distribution to the gY structure function for three 
sets of polarized parton distributions considered in this analysis. 

5 Conclusions 

Charged current events, generated with the PEPSI Monte Carlo program, have been studied 
for a polarized HERA. Detector effects were taken into account. Charged current events allow 
a precise extraction of the structure functions gY in the region of Bjorken-x larger than 0.01. 
Systematic effects of the calorimeter energy scale, parton distributions and detector acceptance 
were studied, with data samples corresponding to a total statistics of 500 pb -1 . The measurable 
asymmetries remain very prominent and survive the detector effects. With simplistic assump- 
tions these structure functions are directly related to the polarized quark distributions Am and 
Ad. Hence the study of asymmetries allows to extract information on these quantities at a Q 2 
scale of 0(1000) GeV 2 . Although present parton density parametrizations do not predict large 
differences, mainly due to the limited information presently available for extracting these densi- 
ties, this measurement is unique and will provide vital data on the flavour decomposition of the 
polarized quark densities. It will be also of interest to check the evolution of these densities at 
low Q 2 as measured by SMC and HERMES |19f to the kinematic region of HERA. 
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